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ABSTRACT 

' We present a method to derive positions of molecular clouds along the lines of sight 

from a quantitative comparison between 2.6 mm CO emission lines and 18 cm OH ab- 
sorption lines, and apply it to the central kiloparsecs of the Milky Way. With some 
simple but justifiable assumptions, we derive a face-on distribution of the CO bright- 
ness and corresponding radial velocity in the Galactic centre without any help of kine- 
matical models. The derived face-on distribution of the gas is elongated and inclined so 
\ that the Galactic-eastern (positive longitude) side is closer to us. The gas distribution 

is dominated by a barlike central condensation, whose apparent size is 500 x 200 pc. 
A ridge feature is seen to stretch from one end of the central condensation, though 
t^J" , its elongated morphology might be artificial. The velocity field shows clear signs of 

' noncircular motion in the central condensation. The 'expanding molecular ring' fea- 

ture corresponds to the peripheral region surrounding the central condensation with 
the Galactic-eastern end being closer to us. These characteristics agree with a picture 
I 1 in which the kinematics of the gas in the central kiloparsec of the Galaxy is under a 

strong influence of a barred potential. The face-on distribution of the in situ pressure 
of the molecular gas is derived from the CO multiline analysis. The derived pressure 
^ ' is found to be highest in the central 100 pc. In this region, the gas is accumulating 

and is forming stars. 
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1 INTRODUCTION 

The behavior of molecular gas, in particular its physical 
conditions and kinematics, in central regions of galaxies 
is key information to understand the star forming activity 
which occurs there. The Galactic centre can be observed in 
much greater detail compared with central regions of other 
galaxies. However, its inevitable edge-on perspective some- 
times complicates the interpretation of the data. In partic- 
ular, a face-on image of the Galactic centre is very hard 
to construct, though such an image would be very help- 
ful to understand its kinematics and to make a compari- 
son with central regions of other galaxies. Attempts have 
been made to construct models of gas kinematics (s ee, e.g., 
iLiszt fc Burtonlll980l : fBinnev et alJll99ll: |Puxlll999l) . Their 
models give more or less reasonable interpretation of some 
aspects of the distribution and motion of the gas, in partic- 
ular its asymmetric position-velocity appearance and some 
outstanding features. Kinematical models can be used to 
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project position-velocity diagrams of molecular /atomic line 
data into a face-on view (see e.g. fCohen fc DentJ Il983l ; 
ISofuelll995l: iNakanishi fc Sofuej I2005T . This is an indirect 
method to investigate the spatial distribution of the gas: 
it would be invaluable if we could derive positions of molec- 
ular clouds independent of kinematical assumptions. The 
molecular content and star forming regions in the Galactic 
centre are strongly confined in the central a few hundred par- 
secs: the region is calle d the central molecular zone (CMZ; 
iMorris fc SerabvrJl996l) . Though the kinematics and face-on 
distribution of the gas within the CMZ are essential infor- 
mation to study the star forming activity in the Galactic 
centre, they are complicated and less well understood than 
the outer region (up to the Galactocentric radius of a few 
kpc) mainly traced by atomic hydrogen line. 

In this paper, we present a method to derive a molecular 
face-on view of the Galactic centre without any help of kine- 
matical models. In § |2] we describe the basic methodology. 
Using that, we draw a face-on distribution of the molecu- 
lar gas from existing data and discuss the resultant face-on 
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view in § |3| Physical conditions of the gas and the validity 
of parameters used in the model are also discussed. 



2 THE METHOD 

The essence of our method lies in comparing emission and 
absorption spectra toward the Galactic centre. We compare 
the CO 2.6 mm emission with the OH 18 cm absorption. Be- 
cause the Galactic centre region itself is an intense diffuse 
18 cm continuum source, strong OH absorption arises pref- 
erentially from the gas that lies in front of the continuum 
regions, rather than gas that lies behind them. On the other 
hand, the CO emission samples the gas both in front and 
back of the continuum sources equally. Thus the OH/CO ra- 
tio carries information on the position of the gas along the 
line of sight relative to the continuum sources. In the follow- 
ing, we devise a method to estimate the position quantita- 
tively by introducing some simple assumptions. 

2.1 Qualitative inspection on molecular gas 
distribution 

We used the Colu mbia-U. Chile 12 CO J = 1 — data 
iBitran et al.lfl997l) and a large-scale absorption survey of 
the OH main l i nes ( 1665 MHz and 1667 MHz) made by 
iBovce fc Cohenl lll994f) . Both data were binned over suc- 
cessive 5kms _1 velocity ranges. The CO data, for which 
the original resolution and grid spacing are respectively 9 
arcmin and 7.5 arcmin, were smoothed to 10-arcmin reso- 
lution and resampled with a 12-arcmin grid to match the 
OH data. Hereafter we mainly used the 1667 MHz line for 
the OH absorption. However, because of the large velocity 
spread of the molecular clouds in the Galactic centre, the 
positive- velocity end of the 1667 MHz line blends together 
with the negative-velocity end of the 1665 MHz line. Thus 
we replaced the ulsr > 100 km s -1 part of the 1667 MHz 
line with the 1665 MHz data scaled by 1.37, which is the 
mean ratio of the absorption depth between the two lines in 
the velocity range 25 ^ vlsr ^ 100 km s -1 . 

Figure shows the processed spectra of the OH ab- 
sorption and the CO emission at (£,b) = (-0?2,0?0). This 
figure demonstrates that the OH/CO ratio varies signifi- 
cantly as a function of the radial velocity, ulsr- For ex- 
ample, we may draw attention to the velocity components 
near «lsr — — 130kms _1 and ulsr — 160kms _1 (shad- 
owed in Fig. 0, both of which belong to the so-called 'ex- 
panding mole cular ring' (EMR: [kaifu. Kato fc Iguchilll972l 
IScovilldlT972T) . The CO intensities at these components are 
almost the same. This suggests that the amounts of molecu- 
lar gas are similar in each component. On the other hand, the 
OH absorption depths are strikingly different. The negative- 
velocity component shows deep absorption, while there is no 
distinct absorption around the positive-velocity component. 
We immediately deduce from this fact that the negative- 
velocity component is located in front of strong continuum 
source surrounding the Galactic centre, while t he positive- 
veloci ty component is behind it. This logic led Kaif u et alJ 
(Il972l) to conclude that this feature is expanding away from 
the Galactic centre. 

Figure [5] shows the longitude-velocity (£-v) diagram 
of the ratio between the OH apparent opacity (r app = 
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Figure 1. Sample spectra of CO emission (thin line) and OH 
absorption (thick line) at (£, b) = (— 0?2, 0?0). Shadowed velocity 
components demonstrate the difference of cloud position along 
the line of sight (see text). 
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Figure 2. The ratio between the OH apparent opacity and the 
CO J = 1 - intensity at b = 0?0. Contours are T MB (CO) = 
1,2,4,7,10, and 14 [K], The ratio is shown in the region where 
Tmb(CO) exceeds 1 [K]. Hatched velocity ranges were excluded 
from the analysis because of contamination by clouds well outside 
the Galactic centre region. 



Tabs/Tcont; where T a b s and T cont are line absorption and 
continuum antenna temperatures, respectively) and the 
CO J = 1 — line intensity. The ratio has a clear trend; 
it is smaller in I < 0° and ?;lsr > Okms , while it is larger 
in £ > 0° and wlsr < Okms -1 and in the feature with 
extremely large ve locity width at £ ~ 3°, which is called 
Bania's Clump 2 jBanialll977l) . The ratio in the velocity 
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range — 60kms 



< 



wlsr < lOkms 1 is affected by the 



foreground gas in the Galactic disc and does not show the 
intrinsic value of the Galactic centre. Thus the hatched ve- 
locity ranges in Fig.|5|are excluded from the following anal- 
ysis. 



2.2 Deriving distances to clouds 

Based upon the logic shown in the previous section, we ex- 
tend it to quantitative estimation of molecular gas distribu- 
tion. In order to determine the distances to molecular clouds 
quantitatively, we adopt the following four assumptions: 

(1) At a given £, emission observed at each velocity bin comes 
from a single position along the line of sight. 

(2) The CO line intensity Too at a given velocity is pro- 
portional to the amount of molecular gas in unit velocity 
width; and the OH opacity toh (not the 'apparent' one 
but the real one) at a given velocity is also proportional 
to the amount of molecular gas in unit velocity width. Con- 
sequently, toh = ZTco where Z is a constant. 

(3) The excitation temperature of OH [T cx (OH)] is uniform. 

(4) The 18 cm continuum emission is optically thin and 
arises from a distributed volume emissivity, j(r) [r is 
the Galactocentric radius], that is axisymmetric about the 
Galactic centre as modelled in § 12.31 

Now Z and T ex (OH) are unknown parameters. How to 
determine them is described in § 12.41 validity of them and 
the assumption (1) are discussed in § 13.31 

When a cloud whose OH opacity is toh is located at 
s — so (s is the position along the line of sight), the cloud 
absorbs the continuum intensity behind it, i.e., the contin- 
uum emissivity integrated along the line of sight behind the 
cloud, J"° j(r)ds using the assumption of optically thin 
continuum emission. The absorption depth is written as 
/[l-exp(-roH)][/!^ j(r)ds-T cx (OH)] where / is the beam 
filling factor of OH absorbing gas. The total continuum in- 
tensity observed by us is J j(r)ds, and the apparent opac- 
ity T app of the cloud is expressed as 



/ [1 - cxp(-toh) 



P j(r)d S -T B (OH) 



Iliads 



(1) 



Since T apP and Too (and thus toh = ZTco) are known 
through the observations, we can obtain the value of 
J_° j(r)ds if we know /, Z, and T ex (OH). Then so is derived 

from the value of j(r)ds using a distribution model of 
j(r). 



2.3 Modelling the distribution of continuum 
emissivity 

We assume that the 18 cm continuum emissivity j(r) around 
the Galactic centre is described as a sum of several axisym- 
metric Gaussians: 
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Figure 3. The schematic relation of geometrical parameters in 
the face-on view of the Galactic centre region. 



Figure|3]shows the schematic relation of geometrical param- 
eters. The Galactic longitude I is in degrees. Here r, s, and a% 
are in units of projected distance corresponding to 1° at the 
distance to the Galactic centre: i.e., 150 pc at 8.5 kpc. The 
observed longitudinal distribution of the continuum bright- 
ness T cout {£) can be written as 



T mnt {l) 



j(r) ds. 



(4) 



Since the continuum emissivity beyond the solar circle 
should be negligible compared with that in the centre, 
j{ r ) ds can be replaced with /_ j(r) ds. Thus 
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exp 
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ds (5) 



(6) 



Now Oj and at can be drawn so that the Eq. ||SJ repro- 
duces the observed longitudinal distribution of the contin- 
uum brightness. It is found that observed continuum distri- 
bution is well fitted by up to three components (N = 3). 
For b = 0?0, we obtain ai = 98.3, a 2 = 34.0, a 3 = 11.4; 
<T\ — 0.120, (72 = 0.677, 03 = 7.18. Parameters ai are in 
Kelvins, and ai are in degrees. Figure [1] shows the result of 
our fit. The model reproduces the observed data quite well. 



Since ai, <T2 are small enough (< 100 pc) and m, ai are 
large, the continuum distribution due to the first (i = 1) 
and the second (i = 2) components is well defined. In some 
;alaxies seen in nearly face-on perspective, such as M83 
Ondrechenlll985l) and IC 342 JCrosthwaite. Turner fc Hoi 
2000), A ~ 20 cm continuum emission in the central kilo- 
parsec consists of strong central source and extended emis- 
sion seen in resolutions of hundreds of parsecs. For such 
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Galactic Longitude [degrees] 

Figure 4. Longitudinal distribution of the 18 c m continuum in 
units of antenna temperature at b = 0? by iBovce &: Coher] 
11994) (open circles) and a fit by 3 Gaussians (solid line). The 
fitting function is T con t = o» exp (— £ 2 /(2<Ti 2 )) where a\ = 

98.3 [K], a 2 = 34.0 [K], a 3 = 11.4 [K]; oi = 0?120, cr 2 = 0?677, 
o~3 = 7? 18. 



continuum distribution, axisymmetric emissivity distribu- 
tion is a good approximation in the first order. On the 
other hand, there are galaxies whose central A ~ 20 cm 
emiss ion is dominated by ring structure , such as NGC 
1097 llHummel. van der Hulst fc Keellll987D and NGC 1365 
iSandavist^5rsate^^^Lln^ladT^95lK ~Tliis is not the case 
for the Milky Way Galaxy because the 18 cm continuum 
brightness in the Galactic centre, our edge-on view, has an 
intense central peak with a monotonic decrease away from 
the centre (see Fig.[IJ. The fact that the central peak, which 
should originate in the Sgr A region, is seen through the 
foreground medium supports our assumption that the 18 cm 
continuum is optically thin. 

On the other hand, because of large 03 and small 0,3, 
the third (i = 3) component should suffer from possible 
non-axisymmetric distribution on the largest scale and/or 
individual sources; furthermore, small dr cont /ds causes large 
errors in the derived positions of clouds. Therefore, in longi- 
tudes where the contribution from the first and the second 
components is negligible (i.e., \£\ > 1?5), the positions of 
clouds obtained using this model are rather uncertain. 



Eland© Trials have been done for Z = 0.04 to 0.70 [K _1 ] 
and T cx (OH) = to 10 [K] \ 

We have chosen an appropriate set of (Z, T cx (OH)) so 
that the following three conditions are satisfied: 

(1) The resultant face-on distribution of the CO brightness 
is not too asymmetric between the near and far sides with 
respect to the centre, 

(2) The features extending above and below the Galactic 
plane are placed in similar face-on positions at different lat- 
itudes, and 

(3) Most of the CO emission has a solution for the position 
s . 

In order to check the condition (1), the CO emission 
within \i\ ^ 1° is treated. The percentage of the CO emission 
located in the near side (so > 0) is displayed in Figure [7Ji. 
The central molecular e mission has a well-k nown asymmetry 
with respect to £ = 0° : iBallv et alj ([1988) mentioned that 
about three-fourths of molecular ( 13 CO and CS) emission 
comes from positive longitudes. Therefore asymmetry to a 
similar degree along the line of sight is acceptable. Sets of 
large Z and small T cx (OH) are rejected since they would put 
most of the emission on the far side. 

A cloud complex around (•£, vlsr) — (1°3, +100 kms - ) 
[hereafter the '1?3 region'] is chosen to check the condition 

(2) . At first we have derived the centres of emission along 
the lines of sight in b = 0?0, ±0?2, ±0?4. Then the stan- 
dard deviation of the centres is calculated: Figure [7J} shows 
the results. The standard deviation becomes small when we 
adopt Z > 0.10 [KT 1 ] and T OX (OH) ~ 4 [K]. The gas at high 
latitudes goes farther with respect to that at b = 0?0 if we 
adopted smaller values for T OX (OH); the reverse is also true. 

For the condition (3), we deal with Clump 2. Because of 
its large OH/CO ratio, Clump 2 is apt to have no solution; 
i.e., the large absorption depth cannot be reproduced even 
if the cloud is located just in front of us. Figure^ shows the 
percentage of the CO emission which can be placed within 
-3750 < s < 3750 pc for the used sets of (Z,T ex (OH)). If 
we use small Z or large T cx (OH), significant emission is lost. 

By combining these conditions, we have chosen Z — 
0.15±0.03 [K _1 ] andT cx (OH) =4±1[K]. As noted in SlO 
the derived distances to the clouds are uncertain in \£\ > 1?5. 
Thus the constraint on the parameters given by condition 

(3) using Clump 2 (£ ~ 3°) is weaker than those given by 
conditions (1) and (2). However, parameter space allowed by 
condition (2) is involved in that given by condition (3). Thus 
we can obtain (Z,T cx (OH)) = (0.15 ± 0.03 [K" 1 ], 4 ± 1 [K]) 
using the conditions (1) and (2) alone; condition (3) is then 
automatically satisfied. 



2.4 Choice of the Z and T GX (OH) values 

Following the procedure, we can draw a face-on (x-y) dis- 
tribution of the molecular gas by putting each data point of 
the £-v diagram (Fig. |2J on to the x-y plane with a projec- 
tion of (£,s) — > (x, y). However, /, Z, and T cx (OH) are still 
unknown. Here we assume / = 1: the validity of this assump- 
tion is discussed in § 13.3.11 There is no bottom-up scheme 
to determine Z and T OX (OH). To determine Z and T cx (OH), 
we employ a trial-and-error scheme, making face-on maps 
at b = +0?4, +0?2, 0?0, -0?2, and -0?4 with various values 
of Z and T cx (OH). The data used, i.e., the OH/CO t-v dia- 
grams and the continuum distribution, are shown in Figures 



3 RESULTS AND DISCUSSION 

3.1 Distribution and kinematics of the gas 

3.1.1 Overall structure 

Figure |8] shows the resultant molecular face-on view of the 
Galactic centre region at b = 0?0 seen from the direction 
of the north Galactic pole. The CO brightness of each pixel 

1 T ex (OH) is measured as an excess over the cosmic microwave 
background (CMB). 
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Figure 5. Longitude- velocity distribution of the OH/CO ratio, shown in the same way as Fig, l2l 



in Fig. |21 is projected on to the x-y plane and smoothed by 
a = 20 pc Gaussian. Figures [HJi and [HJ> show the distribu- 
tions of CO brightness and corresponding radial velocity, 
respectively. Dashed lines in the Figure indicate t = ±1?5: 
as mentioned in § 12.31 the obtained face-on maps are less 
reliable outside them. 

The resultant CO distribution mainly consists of a cen- 
tral condensation (corresponding to the CMZ) and a thin 
ridge feature at (x, y) ~ (400 pc, —1000 pc) that stretches 
almost along the line of sight from the Galactic-eastern (pos- 
itive x in Fig. |SJ end of the central condensation. As seen in 
Fig- HI high OH/CO ratio is more widespread in the positive 
longitudes. This comes about because the face-on CO dis- 
tribution is inclined so that the gas in the positive x (corre- 
sponding t o positive I) l i es clo ser to us, as qualitatively sug- 
gested bv lCohen fc Fewl (£l976). It has long been argued that 
the Milky Way is a barred galaxy. Recent studies agree with 
a picture in which the Galactic-eastern ( positive longitude) 
end of the bar is closer to us (see, e.g.. iMorris fc Serabvnl 



Il996t iGerhardl Il998t and references therein). Our face-on 
map shows, as a whole, the same trend as this picture. 

Some barred galaxies, su c h as NGC 1097 
JOndrechen fc van der Hulsti Il983t iHummel et~ai] Il987fl 
and M83 lbndrechenlll985l) . show A ~ 20 cm continuum 
enhancements in their bars; in particular, in the dust 
lanes. If the Milky Way is barred and the 18 cm continuum 
emissivity is elongated along its bar like these galaxies, 
then the actual distribution of molecular gas would be more 
inclined than that obtained here. Therefore the trend of 
inclination is very certain. 



3.1.2 Central condensation 

The central condensation dominates the CO emission in the 
Galactic centre region. It is elongated: its apparent size is ap- 
proximately 500 x 200 pc, and should be comp ared with 'twin 
peak s' in central regions of barred galaxies llKennev et alJ 
1992). The minor axis length of the condensation might be 
smaller since the face-on map involves positional errors along 
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Figure 6. Longitudinal distribution of the 18 cm continuum brightness and model fits shown in the same way as Fig.|4] The fit parameters 
a; and ai (i = 1, 2, 3) are in Kelvins and degrees, respectively. 
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Figure 7. (a) Percentage of the CO emission within \£\ 1° which is located in the near side. The unit is per cent, (b) The standard 
deviation of the centres of CO emission in a cloud complex (1?0 ^ I $C 1?6, 40 i>LSR Ss 130 km s ) at b = 0?0, ±0?2, ±0?4 in parsecs. 
(c) Percentage of the CO emission in Clump 2 (2? 6 < t < 3? 4, 0?0 < b < 0?4, 30 < v LSR < 200 kms" 1 ) which can be placed within the 
inner 3750 pc. The unit is per cent. Desirable parameter space is shadowed. 



the lines of sight as noted later in § 13.1.31 The major axis of 
the condensation is inclined with respect to the line of sight 
(x = 0) by ~ 70° so that the Galactic-eastern side is closer to 
us. This angle does not change significantly in the acceptable 
parameter space of (Z, T cx (OH)). The condensation includes 
Sgr A (Z ~ 0?0), Sgr B [l ~ 0?6), and Sgr C {I ~ -0?5) 
molecular cloud complexes and the 1?3 region. The resul- 
tant face-on distribution of radial velocity (Fig. |Hp) clearly 



shows that the gas motion in the condensation is strongly 
noncircular. The gas in the far-side has larger receding ve- 
locity, while the gas in the near-side is rather approaching. 
This velocity structure can be explained if the gas orbit is 
elongated along the major axis of the central condensation. 
It is a clear sign of a bar: a similar trend is often seen in both 
numerical simulations of gas kinematics in a b arred potential 
and observations of barred galaxies (see, e.g.. lAthanassoulal 
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Figure 8. The resultant molecular face-on view of the Galactic centre at b = 0?0. The parameter set (Z,T ex (OH)) = (0.15K" 1 ,4K) 
is used. Distribution of CO J = 1 — emission (normalized by the peak value) is shown in (a). Contours are 0.01, 0.02, 0.04, 0.08, 
0.16, 0.32, and 0.64 of the peak value. Corresponding radial velocity is shown in (b). Distribution of the velocity is shown in the region 
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Thick contour is i>lsr = Operas -1 ]. The solar system is located at (x,y) = (0, —8500). White crosses in both panels show the centre, 
(x,y) = (0,0). Dashed lines indicate I = ±1?5. 



Il992l: iLindblad. Lindblad fc Athanassoulalll99d) . This fact 
supports the arguments that the Milky Way is barred. 



3.1.3 Bania's clump 2 

The ridge, which lies at (x,y) — (400 pc, — 1000 pc), cor- 
responds to the feature called Bania's Clump 2. Clump 2 
has long been noticed because of its peculiar characteristics: 
widespread velocity s tructure and large latitudinal extent. 
IStark fc Banial il986t) considered that Clump 2 can be un- 
derstood as an end-on projection of a dust lane or an inner 
spiral arm. iFuxl lll999l) argued that, however, a dust lane 
which associates with the Galactic-eastern side of the bar 
does not corresp ond to the Clump 2, b ut to the 'connect- 
ing arm ' fe ature (iRoueoor fc Oortll960D based on numerical 
simulation. IFuxl lll999l) discussed Clump 2 and interpreted 
it as gas clouds which are crossing a dust lane shock. 

In our face-on map, the ridge corresponding to Clump 
2 resembles offset ridges (or dust lanes) often seen in barred 
spiral galaxies. Naively speaking, the entire structure, a cen- 
tral condensation and a possible offset ridge-like feature, is 
similar to those seen in centr al regions of barred galaxies 
(see, e.g., iKennev et al] ll992!). It also mimics numerically 
simulated distrib utions of materials o rbiting in a barred po- 
tential (see, e.g.. lAthanassoulalll992l) . 



However, the elongation of the ridge along the line of 
sight might be artificial: the face-on maps obtained here in- 
evitably involve positional errors along the lines of sight. 
The nominal error reduces around the centre since the gradi- 
ent of continuum brightness (dT cont /ds, namely, continuum 
emissivity j) becomes steeper; and vice versa. The la devia- 
tions in Too and T app typically cause ~ 200 pc error around 
Clump 2, while the error in the central condensation is a 
few tens of parsecs. A localized source of radio continuum 
emission embedded in the molecular cloud could also lead, 
in our model, to an apparently more widespread distribution 
of the cloud in the derived face-on map than the actual one. 

Therefor e our results cannot de term ine which o f the in- 
terpretations |Stark fc Banial dl986l) and lFuxl 119991) ] is rele- 
vant. Nevertheless, it is certain that Clump 2 lies in front of 
the centre, since the OH/CO ratio there is remarkably high. 



3.1.4 'Expanding Molecular Ring' feature 

The so-called EMR feature was considere d as an expanding 
(and rotating ) ring in the early works llKaifu et al.l 11973 : 
IScovilldll972l) . as obvious from its name, based on its tilted 
oval t-v appearance. Later, another interpretation that the 
apparent expa nsion results from e longated orbit of the gas 
was proposed. iBinnev et al.l il99lt) claimed that the l-v lo- 
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Figure 9. Longitudinal distribution of CO J = 1 — integrated 
intensity (thin line) and OH integrated apparent opacity (thick 
line) for negative- velocity component of the EMR. The velocity 
range for the integration is —300 ^ i>LSR [kms - 1 ] $C —65 + 5£ 
(I is in degrees): approaching toward us faster than the hatched 
velocity ranges in Fig. l2l 



cus of the EMR, a parallelogram rather than an ellipse, 
can be understood as a projection of the so-called x\ or- 
bit (|ContoDoulos fc Mertzanideslll977T) in a barred poten- 
tial. As noted in the present (see § 12. II and past (see, e.g., 
iKaifu et al]ll972l) works, the negative-velocity component 
of the EMR lies in front of the centre, while the positive- 
velocity component is located behind of the centre. Both 
interpretations, an expanding ring and an orbit in a barred 
potential, were constructed to match this geometry. 

In the face-on map obtained in this work, the EMR (in 
particular negative-velocity component) is not clearly sep- 
arated from the central condensation. We investigate the 
origin of the EMR, an energetic expansion or an orbital 
crowding, based on the OH/CO ratio. We would pay atten- 
tion to the negative-velocity component of the EMR, which 
is prominently seen in l-v diagrams at b = 0?0, — 0?2, — 0?4 
(Fig. [SJ . Higher ratio is more extensively seen in the pos- 
itive longitudes. For example, we present the longitudinal 
distribution of the CO emission and the OH absorption at 
b — — 0?2 in Figure [5] The Figure shows the trend that 
the OH absorption is deeper in the positive longitudes com- 
pared with the CO emission. This suggests that the EMR, at 
least its negative-velocity component, is inclined so that the 
positive-longitude end is c loser to us. This geo metry agrees 
with the interpretation bv lBinnev et alJ (Il99ll) rather than 
the axisymmetric 'expanding ring'. 



3.2 Physical conditions 

ISawada et alJ (2001) have investigated £-v distribution of 
physical conditions of molecular gas in the Galactic centre 
using the CO intensity ratios 12 CO J = 2-l/ 12 CO J = 1-0 
[i? 2 -i/i-o( 12 CO)] and 13 CO J = 2 - l/ 12 CO J = 2-1 
\R\zi\i{J = 2—1)]. In this subsection we discuss the face-on 
distribution of CO intensity ratios and physical conditions 
of the gas derived from the ratios. 

In the face-on map derived in § 13.11 (Fig. ISt. the Sgr B 
cloud complex at (•£, dlsr) — (0?6, 50 kms -1 ) is mapped on 



to (x,y) ~ (90pc, — lOOpc). That is highly offset from the 
other part of the central condensation toward us, and the 
isovelocity contours are unnaturally curved around there. 
This is caused by the intense discrete continuum sources 
Sgr Bl and B2. Therefore we reassigned the cloud positions 
toward £ = 0?6 by interpolating the velocity field from £ — 
0?4 and 0?8: i.e., data at each velocity bin toward £ — 0°6 
are placed on to the position of the corresponding velocity in 
the interpolated velocity field. The resultant face-on maps 
of CO and corresponding radial velocity are presented in 
Figures [TUb and llOb . respectively. 

The velocity field in Fig. I10t > can be used to con- 
vert an £-v plot of any parameter (intensity of other lines, 
line intensity ratios, etc.) to a face-on view. We pro- 
jected the 12 CO J = 2 - 1 and 13 CO J = 2 — 1 data 
JSawada et alJl200ll) at b = 0?0 on to the x-y plane. Fig- 
ures 110b and HOb i show the face-on distribution of inten- 
sity ratios i?2-i/r-o( 12 CO) and R\%/\2{J = 2 — 1), respec- 
tively. On one hand, i? 2 -r/i-o( 12 CO) is high (1.0-1.2) al- 
most all over the central condensation. On the other hand, 
high Rn/i2{J = 2 — 1) (i.e., 0.10-0.15) is seen mainly within 
a radius of ~ 100 pc around the centre. 

Using one-zone larg e velocity gradient (LVG: see, e.g., 
iGoldreich fc Kwanlll974l) analysis, ISawada et all l200ll) de- 
rived the physical conditions of the gas from these intensity 
ratios. Though the derived physical conditions depend on 
the assumed gas kinetic temperature Tk, the thermal pres- 
sure p/k = n(H2)Tk ["-(H2) is the derived number density of 
molecular hydrogen] was found to be almost independent of 
the assumed T^. 

Figure fTTl a shows the distribution of pressure n(H2)Ik 
obtained in the same way. The highest pressure, 10 4 ' 9 - 
10 5 ' 1 cm -3 K, is seen within the central 100 pc. It is note- 
worthy that the highest pressure is found roughly symmet- 
rically around the centre, whereas the molecular mass dis- 
tribution is asymmetric in the way that the larger fraction 
of the CO emission comes from the positive longitude (see 
Fig. EOt). The central region corresponds to the 'high pres- 
sure r e gion' in the l-v diagram des c ribed in ISawada et alJ 
j200ll) . lLaunhardt. Zvlka fc Mezgerl (120021) investigated the 
inner hundreds of parsecs of the Galaxy using IRAS and 
COBE data, and deduced that gas and dust within the galac- 
tocentric radius of 120 pc (referred to as 'inner warm disk' 
in their paper) are warmer than those at larger radii ('outer 
cold torus'). The 'high pressure region' fou nd in our analysis 
may c orrespond to the 'inner warm disk' in lLaunhardt et al] 
( 2002). The 1?3 region shows slightly lower pressure than the 
average of the central condensation; and the EMR has even 
lower pressure. 

A schematic illustration of th e majo r molecular fea- 
tures is shown in Figure ITTb . Sofuc (1995) identified a pair 
of arm-like features in the cen t ral co ndensation from the 
13 CO data taken by Ball y et alJ (119871) . though these 'arms' 
are not clearly separated in our face-on map. It may be 
due to insufficient spatial resolution of the present anal- 
ysis and to deviation from the assumed smooth, axisym- 
metric continuum emissivity because of embedded discrete 
radio continuum sources. These arms and the high pres- 
sure region occupy a similar location in the l-v diagram 
ISawada et alJl200lD . loka et al.l <1996fl argued that forma- 
tion of massive stars may have been taking place contin- 
uously or intermittently in this region for more than 10 8 
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Figure 10. Face-on distribution (at b = 0?0) of (a) 12 CO J = 1 — 0; (6) Radial velocity; (c) 12 CO J = 2 — l/ 12 CO J = 1 — intensity 
ratio; and (d) 13 CO J = 2 — l/ 12 CO J = 2 — 1 intensity ratio. Values are shown in the region where the 12 CO J = 1 — emission is 
more intense than 2 per cent of the peak. Contours show the 12 CO J = 1 — distribution: levels are 0.01, 0.02, 0.04, 0.08, 0.16, 0.32, 
and 0.64 of the peak. 



years based upon comparisons between their CO data and 
the l-v distributions of H109a recombinati on line emission 
JPauls fc Mezgerlll975l) and OH/IR stars (iLindavist et alJ 
Il99ll) . and called it the 'star-forming ring'. We consider that 
concentration of gas may occur to form high pressure (high 
density and/or high tem perature) clouds in this region, from 
which stars form. Mizutani et al.l il994T) observed the [C n] 
158 urn line, which is considered to trace photodissociation 
regions, within \l\ ^ 0?7 and found that the emission has 
peaks at Sgr Bl and Sgr C, in addition to Sgr A region. This 
might suggest that the star-forming activity is enhan ced in 
the Sg r B and C regions. In the two-arm model by ISofud 
( 1995), Sgr B and C are both located on the leads of these 
arms. In his analysis, 'Arm I' (which contains Sgr B) with 
lower (rather approaching) radial velocity is in the near side; 
'Arm IF with Sgr C is in the far side. In our face-on map, 



the far side of the central condensation shows larger reced- 
ing velocity (see S I3.1I . This geometry is consistent with that 
suggested bv lSofuej <ll995l) . therefore Sgr B and C might be 
actually at the leads of these arms. The geometry that the 
sites of active star formation are located at the leads of inner 
arms may suggest a time lag between the arrival of the gas 
into the cent ral orbit and the b eginning of s t ar for mation as 
discussed bv lKohno. Kawabe fc Vila-Vilard (Il999l) for NGC 
6951. 

The 1?3 region lies at the Galacti c-eastern end of the 
central condensation. loka et al . (1998b) noted that the large 
velocity extents and the fluffy structure of the clouds and fil- 
aments in this region may b e a sign that these clouds are 
formed by large-scale shocks. iHiittemeister et al.l (Il998tl re- 
ported that SiO in this region is exceedingly abundant and 
its emission arises from hot, thin gas; which also indicates 
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Figure 11. (a) Face-on distribution of gas thermal pressure derived from one-zone LVG analysis at b = 0?0. Values are shown in the 
region where the 12 CO J = 1 — emission is more intense than 2 per cent of the peak, (b) Schematic illustration of some molecular 
features. Contours in both panels are the same as those in Fig. 1101 
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Figure 12. A picture of molecular gas in the Galactic centre 
region proposed on the basis of the results from this work. Linear 
scales are approximate. 



the existence of shocks. From their arguments and its geom- 
etry, we suppose that the gas in the 1?3 region has recently 
arrived in the central region and collided with the central 
condensation. The extremely large scale height of this re- 
gion compared with other parts of the central condensation 
supports our idea: the gas which has just fallen into the cen- 
tral gravitational field has not relaxed yet. The envelope of 
the 1?3 region toward the Galactic-east and smaller radial 
velocity is placed closer to us than the main body of the 1?3 
region in the face-on map. This may support our idea that 
the 1?3 region has just bumped into the central condensa- 
tion: the envelope may be a t ail of the cloud that has not 
yet reached the condensation. ISakamoto. Baker fc Scovillel 



( 2000) suggested that 'episodic fueling' is occurring in NGC 
5005: the 1?3 region might be a result of such kind of cloud 
infall. Our results are schematically summarized in Figure 

We have found that the major axis of the central con- 
densation is inclined by an angle of ~ 70° with respect to 
the line of sight. Therefore, the position angle between the 
major axes of the large-scale stellar bar and the central 
condensation is ~ 40°-50° (central condensation is lead- 
ing the bar), if our viewing angle of the bar is ~ 20°-30° 
fsee [ Morris fc Serabvnlll99fJ: iGerhardl Il999l iDeeuchi et ail 
2002, and references therein). In some barred galaxies, 
there are central mo lecular gas concentrations ('twin peaks', 
iKennev et"ai] [l992'). In a face-on projection of high resolu- 
tion images, some of these concentrations tend to consist of 
a pair of arms and their major axes are inclined with re- 
spect to the bar major axes b y moderate a ngles: such as 
IC 342 bv llshizuki et al] (Il990|). M101 etc bvlKennev et alJ 
Jl992i). NGC 1530 bvlRevnaud fc Downesl Jl999h, NGC 695 1 
bv iKohno et~al"l Jl999t). NGC 5383 bv ISheth et al] feOOfl . 
and NGC 4303 bv lSchinnerer et al] 12002^ . Thev are similar 
to our face-on view of the Galactic centre. Hydrodynami- 
cal simulations of the barred galaxies also produce similar 
gas distributions: e.g ., simulations for NGC 1365 made by 
iLindblad et al] il996t) produce a central condensation which 
is inclined from the bar a xis with a similar a ngle as that 
we have obtained, 40°-50°. lBinnev et al] il99ll) interpreted 
the longitude-velocity feature that corresponds to the cen- 
tral condensation as a projection of gaseous materials in X2 
orbits. However the X2 orbits are aligned perpendicular to 
the potential (Fig. 3 in lBinnev et al.lll99l]) . As a result, the 
sense of inclination with respect to the line of sight are op- 
posite between our and their interpretation. This difference 
may be due to the fact that the xi or bits perpend icular to 
the bar are stellar (collisionless) orbits. Wada| dl994T) showed 
from his damped-orbit model that gaseous orbits do not 
align with respect to the bar but their major axes gradually 



Molecular Face- on View of the Galactic Centre 11 



1500 



1000 



500 



& 



-500 



-1000 



-1500 



Figure 13. Contours of integrated continuum emissivity 
j(r)ds (see Eq. 0) at b = 0?0 based upon the 3 Gaus- 
sians model. Contour levels are 2.5, 5.0, 7.5, 10.0, and 20.0 K 
in units of brightness temperature. A molecular cloud located on 
this plane should be seen in emission if T ex (OH) is higher than 
the integrated continuum emissivity shown by the contours. Co- 
ordinates are the same as Fig. [S] The centre (x,y) = (0,0) is 
shown by a cross. Dashed lines indicate I = ±1?5. 
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Fieure tTol shows contours of integrated continuum emis- 
sivity along lines of sight, J*_^° j(r)ds in Eq. at b — 0?0. 
It is expected that the gas behind the contour which is equal 
to T ex (OH) is seen in emission, while the gas in front of the 
contour is seen in absorption. However, no significant OH 
emission is observed in this region. Thus, if T cx (OH) > 6 [K], 
the face-on gas distribution is almost restricted to near side 
of the Galactic cent re, y < 0. We consider it is unrealistic. 
On the other hand, iBovce fc Cohenl (I1994T) noted that the 
OH absorption was only detected at positions where the con- 
tinuum temperature exceeded 5 K of antenna temperature 
(7.5 K of brightness temperature). Lower limit of excitation 
temperature of OH is provided from this fact: T ox (OH) is 
higher than the integrated continuum emissivity if the OH 
lines are not seen as absorption. For a cloud in the Galactic 
centre region, integrated continuum emissivity at the posi- 
tion of the cloud is typically a half of the observed bright- 
ness temperature. Thus T ex (OH) should not be far below 
(7.5/2) K. 

Again using Eq. Q , we can show that the upper limit of 
T cx (OH) is described as (l-r app ) /*® j(r)ds: T cx (OH) gets 
the highest if the cloud is located in the solar neighborhood. 
Toward Clump 2, this value goes down to 7K, which gives 
a similar constraint to T ox (OH) < 6 [K] derived above. On 
the other hand, our trials have shown that almost all clouds 
are placed on the far side of the Galactic centre (y > 0) at 
high Galactic latitudes if T ex (OH) is lower (< 3 [K]). 

In conclusion, 4 K is a proper estimate of the excitation 
temperature. 



changes with the Galactocentric radius. He showed that the 
orbits at the Inner Lindblad Resonances lead the bar by 45° : 
the trend agrees with the gas distribution we have obtained. 



3.3 Validity of parameters 

3. 3. 1 Beam filling factor 

We have assumed that /, the beam filling factor of the OH 
absorbing gas, is equal to unity. From Eq. Q, we can show 
/ > Ta PP : / gets the lowest if the cloud is located in the 
solar neighborhood and has infinite opacity. For some lines 
of sight, the apparent opacity r app , which is the lower limit 
of /, is rather large: 0.64 toward Sgr B2; 0.35 toward the 
1?3 region and Clump 2. This fact supports that the beam 
filling factor is large and can be reasonably replaced with 

unity. 

ISawada et alJ (1200 if) estimated the beam filling factor 
of CO J = 1 — emission to be 0.4-0.7 based upo n the LVG 
analyse s and high resolution CO data taken by lOka et all 
(1998b). Moreover, since it is expected that the OH absorp- 
tion arises also from less dense cloud envelope compared 
with the CO emission, the beam filling factor of OH absorb- 
ing gas is at least similar and can be even larger. This again 
suggests that / is nearly unity. 



3.3.2 Excitation temperature 

We have assumed that the excitation temperature of OH, 
T cx (OH), is uniformly 4K above the C MB. It is rep orted 
that T OX (OH) typically ranges in 5-10 K <Elitzurll992T) . Our 
value, 4K + Tcmb, is close to this. 



3.3.3 The 'Z' factor 

We have derived Z, the factor to convert Tco into toh, in 
§ I2.2l empiricallv. Here we attempt an independent estimate 
of Z. 

There are three previously-known relations. (1) It is 
considered that the column density of molecular hydro- 
gen JV(H2) is proportional to the CO J = 1 — inte- 
grated intensity. The conversion factor (the so-called X- 
factor) is measured to be about 2 x 10 20 [cm -2 (Kkms -1 ) -1 ] 
for molecular clouds in the Ga lactic disc (see, e.g., 
iDame. Hartmann fc ThaddeusH200lT) . For the Galactic cen- 
tre clouds, however, it is reported that the X-factor is 
smaller than t hat for the Galac tic disc: i.e., X = several x 
10 19 (see, e.g. JOka et~aflll998all (2) Relative abundance of 
OH to molecular hydrogen, [OHJ/rH^ i s typically several x 
10 -7 (see, e.g.. iHerbst fc Leund Il98gh . (3) The OH col- 
umn density is derived from the OH 1667 MHz opacity r 
as AT (OH) = 2 x 10 14 T OX J rdv. 

Using these relations, the Z factor is written as 



Z\K 
5 x 10 



_r _ [OH]/[H 2 ]X[cm- 2 (Kkms- 



2 x 10 14 
-3 [OH]/[H 2 ] 



T cx [K] 
X 10K 



10- 



10 20 cm" 2 (Kkms" 1 )" 1 T CJ 



(7) 
(8) 



This value becomes several x 10 , which is several times 
smaller than that we have adopted. However Z ~ several x 
10~ 2 does not seem to be realistic, because such a small Z 
cannot reproduce the observed absorption depth for a signif- 
icant fraction of the clouds. This discrepancy can be solved 
if we consider larger X (several x 10 20 ) or higher OH abun- 
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dance (several x 10 6 ). On one hand, the large X-factor is 
unrealistic because the small X-factor in the Galactic cen- 



tre region is consistently suggested in various ways (set 


\- e.g., 


Blitz e^ilj 


1985 


;ICox & Laureiia 


1989; Sodroski et al. 


1995; 


lOka et al.l 


1998e 


On the other hand, it is reported that 



the OH abundance is enhanced behind interstellar shocks 
(see, e.g.. lWardlelll999h . In a high resolution CO map of the 
Galactic cent re, innumerable shells/arcs are seen llOka et alJ 
Il998bl.l200iallbl) . If the shells/arcs are formed by a number of 
explosive phenomena such as supernova, the Galactic centre 
is filled with shocks. In such condition, the OH abundance 
may be higher than that in usual condition. High ioniza- 
tion rate caused by X-ra y also enhances the OH abundance 
jLepp fe DargarnoHl996l) . 

3.3.4 One-to-one correspondence between position and 
velocity 

In § 12.21 we have assumed that a value in the radial velocity 
corresponds to a single position at a given line of sight. Of 
course this assumption is nominally wrong: there are two 
(or more) points in a line of sight with a given radial ve- 
locity. Thus, if the Galactic centre region is almost filled 
with molecular gas, the assumption becomes fairly wrong. 
Molecular gas distribution is, however, rather lo calized. In 
fact, o ur (Fig. |5J and higher-resolution (see, e.g.. lOka et alJ 
1998b) molecular line data show that velocity overlapping 
between major features (Sofue's arms, EMR, etc.) occurs in 
very restricted regions. Therefore the assumption might not 
introduce serious errors into the resultant face-on maps. 

Even if overlapping of clouds is occurring, the results 
do not change very much: let we assume that two individual 
clouds (positions si and S2; si < S2) with the same velocity 
is seen as a single cloud. Position of this 'cloud' is derived as 
so in our method. Using Eq. Q, we can show si < so < S2 
if the H2 column density can be fully traced by Too ■ Thus 
the overall CO distribution derived is not heavily affected in 
this situation. 



4 CONCLUSIONS 

We have developed a method to derive positions of molecu- 
lar clouds along the lines of sight. The method is completely 
independent of any kinematic model and based on observ- 
able data alone; the CO emission line, the OH absorption 
line, and 18 cm continuum distribution. It is applied to the 
central region of the Milky Way to obtain a molecular face- 
on map. The obtained spatial and kinematical structure of 
molecular gas shows the following characteristics. 

— The CO distribution mainly consists of a central con- 
densation and a ridge. 

— The central condensation is elongated and its major 
axis is inclined with respect to the line of sight by ~ 70° so 
that the Galactic-eastern end is closer to us. The gas within 
it shows highly noncircular motion: the gas in the far side 
is receding whereas the gas in the near side is approaching. 
This noncircularity of the gas motion is most likely induced 
by a barred potential. 

— The ridge, corresponding the Bania's Clump 2, lies 
closer to us than the central condensation by ~ 1 kpc. This 



feature apparently mimics offset ridges often seen in barred 
galaxies, although its stretch along our line of sight might 
be artificial. 

— The so-called 'Expanding Molecular Ring' feature does 
not appear as a coherent structure on the resultant face- 
on map. However its negative-velocity side, which lies in 
the near side of the Galactic centre, is inclined so that the 
Galactic-eastern s ide is closer to us. T his is consistent with 
the arguments bv iBinnev et al.l (ll99lT) . that the EMR is a 
projection of X\ orbit in a barred potential, rather than the 
original picture of a ring that is expanding and rotating. 

These results give a new evidence for the existence of 
a bar in the Milky Way Galaxy based on direct distance 
derivation independent of kinematic models. A face-on map 
of thermal pressure derived from comparison of three CO 
transition lines and an LVG analysis is also presented. The 
pressure is distinctly higher (i.e., > 10 5 cm -3 K) within the 
Galactocentric radius of 100 pc, compared with the outer re- 
gion. This high p ressur e region coexists with 'G alactic cen- 
tre a rms' (Sof uel 119951) and 'star- forming ring' ( Ok a et al] 
1996). Concentration of clouds would be occurring due to 
some kind of gas orbit crowding. 



ACKNOWLEDGMENTS 

We acknowledge Leonardo Bronfman for providing us the 
CO J = 1 — data in a computer readable form. This work 
was supported by a Grant-in-Aid for Scientific Research of 
the Ministry of Education, Culture, Sports, Science, and 
Technology 08404009 and 10147202. 



REFERENCES 

Athanassoula E., 1992, MNRAS, 259, 345 
Bally J., Stark A. A., Wilson R. W., Henkel C, 1987, ApJS, 
65, 13 

Bally J., Stark A. A., Wilson R. W., Henkel C, 1988, ApJ, 
324, 223 

Bania T. M., 1977, ApJ, 216, 381 

Binney J., Gerhard O. E., Stark A. A., Bally J., Uchida K. 

I., 1991, MNRAS, 252, 210 
Bitran M., Alvarez H., Bronfman L., May J., Thaddeus P., 

1997, A&AS, 125, 99 
Blitz L., Bloemen J. B. G. M., Hermsen W., Bania T. M., 

1985, A&A, 143, 267 
Boyce P., Cohen R. J., 1994, A&AS, 107, 563 
Cohen R. J., Dent W. R. F., 1983, in Burton W. B., Is- 
rael, F. P., eds., Surveys of the Southern Galaxy. Reidel, 

Dordrecht, p. 159 
Cohen R. J., Few R. W., 1976, MNRAS, 176, 495 
Contopoulos G., Mertzanides C, 1977, A&A, 61, 477 
Cox P., Laureijs R., 1989, in Morris M., ed., Proc. IAU 

Symp. 136, The Center of the Galaxy. Kluwer, Dordrecht, 

p. 121 

Crosthwaite L. P., Turner J. L., Ho T. P. T., 2000, AJ, 119, 
1720 

Dame T. M., Hartmann D., Thaddeus P., 2001, ApJ, 547, 
792 

Deguchi S., Fujii T., Nakashima J., Wood P. R., 2002, 
PASJ, 54, 719 



Molecular Face- on View of the Galactic Centre 13 



Elitzur M., 1992, Astronomical Masers. Kluwer, Dordrecht 
Fux R., 1999, A&A, 345, 787 

Gerhard O. E., 1999, in Merritt D., Valluri M., Sellwood J., 
eds., ASP Conf. Ser. Vol. 182, Galaxy Dynamics. Astron. 
Soc. Pac, San Francisco, p. 307 
Goldreich P., Kwan J., 1974, ApJ, 189, 441 
Herbst E., Leung C. M., 1989, ApJS, 69, 271 
Hummel E., van der Hulst J. M., Keel W. C, 1987, A&A, 
172, 32 

Hiittemeister S., Dahmen G., Mauersberger R., Henkel C, 
Wilson T. L., Martin-Pintado J., 1998, A&A, 334, 646 

Ishizuki S., Kawabe R., Ishiguro M., Okumura S. K., 
Morita K.-L, Chikada Y., Kasuga T., 1990, Nat, 344, 224 

Kaifu N., Kato T., Iguchi T., 1972, Nature Phys. Sci., 238, 
105 

Kenney J. D. P., Wilson C. D., Scoville N. Z., Devereux N. 

A., Young J. S., 1992, ApJ, 395, L79 
Kohno K., Kawabe R., Vila-Vilaro B., 1999, ApJ, 511, 157 
Launhardt R., Zylka R., Mezger P. G., 2002, A&A, 384, 

112 

Lepp S., Dalgarno A., 1996, A&A, 306, L21 
Lindblad P. A. B., Lindblad P. O., Athanassoula E., 1996, 
A&A, 313, 65 

Lindqvist M., Winnberg A., Habing H. J., Matthews H. E., 

1991, A&AS, 92, 43 
Liszt H. S., Burton W. B., 1980, ApJ, 236, 779 
Mizutani K. et al., 1994, ApJS, 91, 613 
Morris M., Serabyn E., 1996, ARA&A, 34, 645 
Nakanishi H., Sofue Y., 2003, PASJ, 55, 191 
Oka T., Hasegawa T., Handa T., Hayashi M., Sakamoto S., 

1996, ApJ, 460, 334 
Oka T., Hasegawa T., Hayashi M., Handa T., Sakamoto S., 

1998a, ApJ, 493, 730 
Oka T., Hasegawa T., Sato F., Tsuboi M., Miyazaki A., 

1998b, ApJS, 118, 455 
Oka T., Hasegawa T., Sato F., Tsuboi M., Miyazaki A., 

2001a, PASJ, 53, 779 
Oka T., Hasegawa T., Sato F., Tsuboi M., Miyazaki A., 

2001b, PASJ, 53, 787 
Ondrechen M. P., 1985, AJ, 90, 1474 

Ondrechen M. P., van der Hulst J. M., 1983, ApJ, 269, L47 
Pauls T., Mezger P. G., 1975, A&A, 44, 259 
Reynaud D., Downes D., 1999, A&A, 347, 37 
Rougoor G. W., Oort J. H., 1960, Proc. Nat. Acad. Sci., 
46, 1 

Sakamoto K., Baker A. J., Scoville N. Z., 2000, ApJ, 533, 
149 

Sandqvist Aa., Jorsater S., Lindblad P. O., 1995, A&A, 
295, 585 

Sawada T. et al., 2001, ApJS, 136, 189 

Schinnerer E., Maciejewski W., Scoville N., Moustakas L. 

A., 2002, ApJ, 575, 826 
Scoville N. Z., 1972, ApJ, 175, L127 

Sheth K., Regan M. W., Vogel S. N., Teuben P. J., 2000, 

ApJ, 532, 221 
Sodroski T. J. et al., 1995, ApJ, 452, 262 
Sofue Y., 1995, PASJ, 47, 527 
Stark A. A., Bania T. M., 1986, ApJ, 306, L17 
Wada K., 1994, PASJ, 46, 165 
Wardle M., 1999, ApJ, 525, L101 



